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ERGODIC PROPERTIES OF VISIBLE LATTICE POINTS 


MICHAEL BAAKE AND CHRISTIAN HUCK 


Abstract. Recently, the dynamical and spectral properties of square-free integers, visible 
lattice points and various generalisations have received increased attention. One reason is the 
connection of one-dimensional examples such as ,^-free numbers with Sarnak’s conjecture 
on the ‘randomness’ of the Mobius function, another the explicit computability of correla¬ 
tion functions as well as eigenfunctions for these systems together with intrinsic ergodicity 
properties. Here, we summarise some of the results, with focus on spectral and dynamical 
aspects, and expand a little on the implications for mathematical diffraction theory. 


Dedicated to Nikolai P. Dolbilin on the occasion of his 70th birthday 

1. Introduction 

Delone sets are important mathematical descriptions of atomic arrangements, and Meyer 
sets are special cases with a rich spectral structure [anaiiT!. Particularly well-studied are cut 
and project sets or model sets, which underly the structure of perfect quasicrystals [HHlIHOl IH], 
It is fair to say that such structures are rather well understood. This is much less so if one 
keeps uniform discreteness but relaxes relative denseness. In fact, one might expect to leave 
the realm of pure point spectrum, at least as soon as one has entropy [7]. However, there 
are interesting examples such as fc-free numbers or visible lattice points that have positive 
topological entropy but nevertheless pure point dynamical and diffraction spectrum. They are 
examples of weak model sets, and deserve a better understanding. Here, we summarise some 
of the known results and add to the structure of their topological and spectral properties. 

The paper is organised as follows. In Section [2l we use the visible points of as a 
paradigm to formulate the results for this case in a geometrically oriented manner and to 
develop our notation and methods while we go along. Section [3] extends the findings to A:-free 
points of n-dimensional lattices, while Section [J] looks into the setting of ,^^-free systems as 
introduced in [I8l[25]. Finally, in Section [5l we analyse an example from the number field 
generalisation of m in our more geometric setting of diffraction analysis. For convenience 
and better readability, the more technical issues are presented in two appendices. 

2. Visible square lattice points 

Two classic examples for the structure we are after are provided by the square-free integers 
(the elements of Z that are not divisible by any nontrivial square) and the visible points 
of a lattice (the points with coprime coordinates in a lattice basis). Since our focus is on 
higher-dimensional cases, let us take a closer look at the visible points of Z^, 

V = V 22 = \ IJ pZ2 = {xe Z^l gcd(x) = 1}, 

p prime 

where gcd(x) = gcd(xi,X 2 ) for x = (xi,X 2 )- Note that, throughout the text, p will denote 
a prime number. The set is illustrated in Fig. [TJ and can also be found in many textbooks 
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i§: 


Figure 1. A central patch of the visible points V of the square lattice 
Note the invariance of V with respect to GL(2,Z). 

including [2], where it is shown on the cover, and [3l Sec. 10.4]. The following result is 
standard; see [3l Prop. 10.4] and references therein for details. 

Proposition 1. The set V has the following properties. 

(a) The set V is uniformly discrete, but not relatively dense. In particular, V contains 
holes of arbitrary size that repeat lattice-periodically. More precisely, given an inradius 
p > 0, there is a sublattice of 1? depending on p such that a suitable translate of this 
sublattice consists of centres of holes of inradius at least p. 

(b) The group GL(2,Z) acts transitively on V, and one has the partition 1? = UmeNo^^ 
of T? into Gh{2,'L)-invariant sets. 

(c) The difference set is V — V = 7?. 

(d) The natural density of V exists and is given by dens(F) = where C, denotes 

Riemann’s zeta function. □ 
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Note that big holes are rare, but important; see O Rem. 10.6] for some examples. The 
interesting fact is that they do not destroy the existence of patch frequencies, though the 
latter clearly cannot exist uniformly. For the natural pair correlation (or autocorrelation) 
coefficients, 

r]{x) := lim 

R-^oo TtR^ ' ' 

one finds the following result; compare [3l Lemma 10.6], [101 Thm. 2], as well as [311 Thm. 7]. 

Lemma 1. For each x G Z^, the natural autocorrelation coefficient r]{x) of V exists, and is 
given by 



where i = n,(i — 2p Ri 0.3226. In particular, with gcd(O) = 0, this also gives the density 

^(o) = np(i-p"') = i/c(2). □ 

The autocorrelation measure of V, 

7 = X] 

is thus well-defined, and is a translation bounded, positive definite measure on by con¬ 
struction. Its Fourier transform 7 exists by general arguments, compare [in Ch. 1.4] or [31 
Rem. 8.7 and Prop. 8 . 6 ], and leads to the following result; see [31 Thm. 10.5] and [131 Thm. 
3]. 


Theorem 1. The natural diffraction measure 7 of the visible points V of the square lattice 
T? exists. It is a positive pure point measure which is translation bounded and supported on 
the points of with square-free denominator, the Fourier-Bohr spectrum of 7 , so 


7 = Iik)dk, 

fceQ^ 

den(A:) square-free 


where den(A:) := gcd{n G N j n/c G I?}. In particular, 7(0) = (1/C(2))^ = 36/71“^, and when 
0 k £ Q'^ has square-free denominator den(A:), the corresponding intensity is given by 


Hk) 



n 

p|den(fc) 



2 


□ 


The essence of this result is the pure point nature of 7 together with its explicit com¬ 
putability via an intensity formula in the form of a finite product for any given fc G 
with square-free denominator. Fig. [2| illustrates the diffraction measure. Note that 7 has the 
symmetry group Z^ xi GL(2,Z). 

An alternative view is possible by means of the Herglotz-Bochner theorem as follows. The 
autocorrelation measure 7 is positive definite on R^ if and only if the function 7 : — > R is 
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Figure 2 . Diffraction 7 of the visible points of . A point measure at k with 
intensity I{k) is shown as a disk centred at k with area proportional to I{k). 
Shown are the intensities with I{k)/I{0) > 10“® and k G [0,2]^. Its lattice of 
periods is 1?, and 7 turns out to be GL(2, Z)-invariant. 


positive definite on Z^; see [3l Lemma 8.4], The latter property is equivalent to the existence 
of a positive measure ^ on = M^/Z^ ~ [0,1)^ such that 

Vix) = [ deiy), 

JT2 

where the connection to 7 is established by = 7 | [ 0 , 1 ) 2 , so that j = g * 5j2, where 5s '■= 
Exes denotes the Dirac comb of a discrete point set S. The finite positive measure g is 
a spectral measure in the sense of dynamical system theory. It is related to the diffraction 
measure by convolution; for background, we refer to UM and references therein. We shall 
return to the dynamical point of view shortly. 

Let us pause to comment on the history and the development of this problem. The arith¬ 
metic properties of V are classic and can be found in many places, including [11123]. The 
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investigation of spectral aspects was advertised by Schroeder in [M], see also [35], by means of 
a numerical approach via FFT techniques. These results suffered from insufficient resolution 
in the numerical treatment, and seemed to point towards continuous diffraction components, 
perhaps in line with the idea that the distribution of primes is sufficiently ‘random’. 

Ten years later, on the basis of a formal Mobius inversion calculation for the amplitudes, 
Mosseri argued in [22| that the diffraction should be pure point rather than continuous, 
thus contradicting the earlier numerical findings of Schroeder. This was corroborated in [3] 
with further calculations on the diffraction intensities (still without proof), which gave the 
formula of Thm. [T] above. Also, a rather convincing comparison with an optical diffraction 
experiment was shown, which clearly indicated the correctness of the formal calculation. The 
first complete proof, with a detailed convergence result with precise error estimates, appeared 
in [BMP], and was recently improved and extended in |31j . on the basis of number-theoretic 
results due to Mirsky |26l[27|. 

Simultaneously, due to the renewed interest in the square-free integers in connection with 
Sarnak’s conjecture, the dynamical sytems point of view became more important, as is obvious 
from Ha El EU HD. Here, the focus is more on the dynamical spectrum, which is closely 
related to the diffraction measure as indicated above, and explained in detail in 13 Ej. 

To explain this, let us define the (discrete) hull of V as 

Xy = {t+ P|t G Z2}, 

where the closure is taken in the product topology induced on {0,1}®^ by the discrete topology 
on {0,1}. This topology is metric |38l El] and is also called the local topology, because two 
elements of {0,1}^^ are close if they agree on a large ball around the origin. Clearly, Xy is 
then compact, where here and below we simultaneously view subsets of I? as configurations. 
In particular, the empty set is identified with the configuration 0, and 1? with 1 this way. 
Since V contains holes of arbitrary size, the empty set is an element of Xy. 

For a natural number m, let ■ m- ^ 1? jnil? denote the canonical projection x i-A [x]m, 

where [x]m = x+mZi^. For a subset X C we denote by Xm its image under this projection 
map. It should however be born in mind that, for elements x G Z^, their images under the 
above map will always be written as \x\m rather than Xm- This is due to the occasional need 
of regarding m in the last expression as an index. Let us recall the following result from [10] . 

Proposition 2 (Chinese Remainder Theorem). |10l Prop. 2] For pairwise coprime positive 
integers mi,m 2 , ■ ■ ■ ,mr, the natural group homomorphism 

r 

i=l 

is an isomorphism. In particular, for xi,X 2 ,... ,Xr G Z^, the simultaneous solutions t G Z^ 

of 

[ijmi — 1 ^ ^ I') 

comprise precisely one coset of mim 2 • ... • mr'^'^ in I?. □ 

Next, let us come to a characterisation of Xy. Let A denote the set of admissible subsets 
A of Z^, i.e. subsets A C Z^ with the property that, for every prime p, A does not contain a 
full set of representatives modulo pi?. In other words, A is admissible if and only if 

\Ap\ < p^ 


(Z 


2 ) 
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for any prime p. Since ^ E A (otherwise some point of V would be in for some prime 
p, a contradiction) and since A is a Z^-invariant and closed subset of {0,1}^^, it is clear that 
Xy is a subset of A. By m Thm. 2], the other inclusion is also true. This was first shown by 
Herzog and Stewart [22] for visible lattice points and by Sarnak [33] for the analogous case 
of the square-free integers. In fact, similar statements hold true for various generalisations 
discussed below; cf. 1311 Thm. 6] for the case of A:-free lattice points. 

Theorem 2. One has Xy = A. □ 


It follows that Xy is hereditary, i.e. 


VXeA: (TcX^TeA), 

and in particular contains all subsets of V. In other words, V is an interpolating set for Xy 
in the sense of [42] . which means that 

Xy|y ■.= {xnv \ X e Xy} = {o, i}^. 

Given a radius p > 0 and a point t E Z^, the p-patch of H at t is 


{V-t)nBp{0), 


the translation to the origin of the part of V within a distance p of t. We denote by A{p) the 
(finite) set of all p-patches of V, and by X(p) = |^(/o)| the number of distinct p-patches of 
V. For a p-patch P of V, denote by C-p the set of elements of Xy whose p-patch at 0 is P, 
the so-called cylinder set defined by the p-patch V; compare [16]. Note that these cylinder 
sets form a basis of the topology of Xy. 

The patch counting entropy of V is defined as 


hpc{V) 


lim 

p^cxD vrp^ 


Note that this differs from the definition in [3I1IM], where, in view of the binary configuration 
space interpretation, we used the base 2 logarithm. It can be shown by a classic subadditivity 
argument that this limit exists. Since Xy is hereditary, it follows that V has patch counting 
entropy hpc{V) at least dens(H) log(2) = (G/tt^) log(2). In fact, one has more. 


Theorem 3. [311 Thm. 3] One has hpc{V) = (G/tt^) log(2). 


□ 


The natural translational action of the group Z^ on Xy is continuous and (Xy, Z^) thus is a 
topological dynamical system. By construction, (Xy,Z^) is topologically transitive mmm, 
as it is the orbit closure of one of its elements (namely V). Equivalently, for any two non¬ 
empty open subsets U and W of Xy, there is an element t such that 


[/n (iT + t) / 0. 

In accordance with Sarnak’s findings [33] for square-free integers, one has the following result. 

Theorem 4. The topological dynamical system (Xy,Z^) has the following properties. 

(a) (Xy, Z^) is topologically ergodic with positive topological entropy equal to (G/vr^) log(2). 

(b) (Xy,Z^) is proximal, and {0} is the unique T?-minimal subset of Xy. 

(c) (Xy,Z^) has no non-trivial topological Kronecker factor [i.e., minimal equicontinuous 
factor). In particular, (Xy,Z^) has trivial topological point spectrum. 
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(d) (Xy,Z^) has a non-trivial joining with the Kronecker system given by K = (G,Z^), 
where G is the compact Abelian group W_p{'^^)p and 7? acts on G via addition of 
i{x) = ([x]p), i.e. g g -\- i{x), with g & G and x G 7?. In particular, (Xy,Z^) fails 
to be topologically weakly mixing. 

Proof. The topological entropy of the dynamical system (Xy, Z^) is just hpciV), so the asser¬ 
tion follows from Theorem[3l cf. [TJ Thm. 1], The topological ergodicity mm will follow from 
the existence of an ergodic full (non-empty open subsets have positive measure) Z^-invariant 
Borel measure on Xy; see Theorems [5] and EKb) below. 

For part (b), recall from Theorem [2] that the hull contains many more elements than the 
translates of V. Nevertheless, one can derive from Proposition (Ha) that every element of 
Xy contains holes of arbitrary size that repeat lattice-periodically. This follows by standard 
compactness arguments from considering a sequence of the form that converges 

in the local topology, via selecting suitable subsequences. In particular, let X,Y G Xy and 
a radius p be fixed. Let tp + P he positions of holes of inradius p in X. Choose p' large 
enough such that Bpr{0) covers Bp{0) F, where T is a fundamental domain of P. Then, 
any p'-ho\e of Y (which exists) contains a p-hole of X. Hence, for any p > 0 and any two 
elements X,Y G Xy, there is a translation t G 7^ such that 

(X -h t) n Bp{h) = (y -h t) n Hp(o) = 0, 

meaning that both X and Y have the empty p-patch at —t. In terms of the metric d on 
Xy [381I311IM] this means that d{X-\-t, Y -\-t) < Ifp and the proximality of the system follows. 
Similarly, the assertion on the unique Z^-minimal subset {0} follows from the fact that 
any element of Xy contains arbitrarily large holes and thus any non-empty closed subsytem 
contains 0 . 

Since Kronecker systems are distal, the first assertion of part (c) is an immediate conse¬ 
quence of the proximality of {Xy,7'^). This also implies that (Xy,7^) has trivial topological 
point spectrum; see | 23 ] for an alternative argument that the non-zero constant function is 
the only continuous eigenfunction of the translation action. 

For part (d), one can verify that a non-trivial (topological) joining |20j of (Xy,Z^) with 
the Kronecker system K is given by 

W:= U (mxn(ZAX)J. 

xeXv p 

Since the Kronecker system K is minimal and distal, a well-known disjointness theorem by 
Furstenberg m Thm. II.3] implies that (Xy,Z^) fails to be topologically weakly mixing. □ 

Following [ml El], the natural frequency v{V) of a p-patch P of K is defined as 

(1) v{V) := dens {{t G 7^ \ {V - t) n Bp{0) = V}) , 

which can indeed be seen to exist. 

Theorem 5. |3Tl Thms. 1 and 2] Any p-patch V of V occurs with positive frequency, which 
is given by 

•'(p) = y n 

j-c(z 2 nBp(o))\p p 


\{'PUF)p\ \ 
p2 J ■ 


□ 
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The frequency function v from m, regarded as a function on the cylinder sets by setting 
i'{C'p) := is hnitely additive on the cylinder sets with 


r&A(p) 


1 

det(Z2)| 


Since the family of cylinder sets is a (countable) semi-algebra that generates the Borel a- 
algebra on Xy (i.e. the smallest cr-algebra on Xy which contains the open subsets of Xy), 
V extends uniquely to a probability measure on Xy; compare m Prop. 8.2] and references 
given there. Moreover, this probability measure is Z^-invariant by construction wherefore we 
have a measure-theoretic dynamical system (Xy,Z^,y). For part (b) of the following claim, 
note that the Fourier-Bohr spectrum of V is itself a group and compare [U Prop. 17]. 


Theorem 6. The measure-theoretic dynamical system (Xy,Z^,y) has the following proper¬ 
ties. 

(a) The 7?-orbit of V in Xy is v-equidistributed, which means that for any function 
f E C'(Xy), one has 


1 


lim 


^ fiV + x)= [ f{X)dHX). 


a;eZ2nSfl(0) 


In other words, V is u-generic. 

(b) (Xy,Z^,y) is ergodic, deterministic {i.e., it is of zero measure entropy) and has pure 
point dynamical spectrum. The latter is given by the Fourier-Bohr spectrum of the 
autocorrelation 7, as described in Theorem^ 

(c) (Xy,Z^,y) is metrically isomorphic to the Kronecker system = {G,l?,p), where 

G is the compact Abelian group )p, the lattice 7? acts on G as above and p is 

the normalised Haar measure on G. 


Proof. For part (a), it suffices to show this for the characteristic functions of cylinder sets of 
finite patches, as their span is dense in C(Xy). But for such functions, the claim is clear as 
the left hand side is the patch frequency as used for the definition of the measure n. 

For the ergodicity of (Xy,Z^ , r), one has to show that 


1 


lim 

R-^00 vri?^ 


,z{{Cr + x)n Gq) = y{Gv)v{GQ) 


xe.'i?r\BR{o) 


for arbitrary cylinder sets Gp and Cq; compare im Thm. 1.17]. The latter in turn follows 
from a straightforward (but lengthy) calculation using Theorem [5] and the definition of the 
measure v together with the Chinese Remainder Theorem. In fact, for technical reasons, it is 
better to work with a different semi-algebra that also generates the Borel cr-algebra on Xy; 
see Appendix!^ for the details. 

Vanishing measure-theoretical entropy, i.e. 

hmeas{V) = lim ^ ^ -y(C-p) log I^(C''p) = 0 , 

^ r&A{p) 


was shown in |31l Thm. 4], which is in line with the results of [7]. Alternatively, it is an imme¬ 
diate consequence of part (c) above. As a consequence of part (a), the individual diffraction 
measure of V according to Theorem [1] coincides with the diffraction measure of the system 
(Xy,Z^,y) in the sense of [5j. Then, pure point diffraction means pure point dynamical 
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spectrum [5l Thm. 7], and the latter is the group generated by the Fourier-Bohr spectrum; 
compare [5l Thm. 8 ] and [H Prop. 17]. Since the intensity formula of Theorem [1] shows that 
there are no extinctions, the Fourier-Bohr spectrum here is itself a group, which completes 
part (b). 

It is well known that = (G, Z^,//) has the same pure point dynamical spectrum as 
(Xy,Z^,i/); compare [13] for the details in the case of the square-free integers. In particular, 
the subgroup of given by the points of H [ 0 , 1 )^ with square-free denominator is easily 
seen to be isomorphic to the direct sum 0pZ^/pZ^, wherefore it is the Pontryagin dual of 
the direct product G = ]]][pZ^/pZ^; cf. |32l Sec. 2.2]. By a theorem of von Neumann [30] . 
two ergodic measure-preserving transformations with pure point dynamical spectrum are 
isomomorphic if and only if they have the same dynamical spectrum. This proves part (c), 
which is a particular instance of the Halmos-von Neumann theorem; cf. [21j . 

Alternatively, the Kronecker system can be read off from the model set description, which 
also provides the compact Abelian group. The general formalism is developed in [ 6 ], though 
the torus parametrisation does not immediately apply. Some extra work is required here to 
establish the precise properties of the measure-theoretic homomorphism onto the compact 
Abelian group. Diagrammatically, the construction looks like this: 


^ Up^p 

u u u 

ML W 

t t t 

X <—)■ (x,t(x)) —)■ t(x) 


Here 

L := {(x, /.(x)) I X G Z^} = {(x, ([x]p)) | x G Z^} 
is the natural (diagonal) embedding of into Z^ x 

lT:=n((Z')p\{[0 U) 

p 

satisfies W = dW and has measure fi{W) = 0^(1 “ = 1/C(2) with respect to the 

normalised Haar measure /U on the compact group Wp{'L^)p- Clearly, one has 

M := M{W) := {x G Z^ | i{x) £ W} = V 

The above diagram is in fact a cut and project scheme: L is a lattice (a discrete co-compact 
subgroup) in Z'^ xW^{Z?)p with one-to-one projection onto the first factor and dense projection 
onto the second factor. This means that H is a weak model set |3]. The corresponding ‘torus’ 
is 

T := {Z^xll{Z\)/L cx ll{Z\, 

p p 

with the Z^-action given by addition of t(x) = ([x]p). A similar construction with translations 
from the group in mind is given in [10] : see also [371 Ch. 5a]. 

The so-called torus parametrisation [28| is the Borel map 

(p: T —>■ Xy, 


given by 


{[yp]p)^M{{[yp]p) + W). 
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Clearly, (p intertwines the Z^-actions. Note that, since = A, one indeed has 

M{{[yp\p) + W) = I?\ [j{yp+pl?) € Xy. 

p 

The map ip fails to be injective. For example, the fibre ip~^{0) over the empty set 0 is easily 
seen to be uncountable. Furthermore, (p is not continuous, since, e.g., ¥^(([(0,0)]p)) = V but 

y 9 ( 1 , 0 ) 0 9 :p(([( 0 , 0)]p ,,..., [( 0 , 0 )]p„_,, [( 1 , 0 )]p„, [( 0 , 0 )] 

Nevertheless, employing the ergodicity of the measure ly, one can show that ip is in fact a 
measure-theoretic isomorphism between the two systems; see Appendix [B] for details. □ 

Let us mention that our approach is complementary to that of m- There, ergodicity and 
pure point spectrum are consequences of determining all eigenfunctions, then concluding via 1 
being a simple eigenvalue and via the basis property of the eigenfunctions. Here, we establish 
the y-genericity of V and the ergodicity of the measure i' and afterwards use the equivalence 
theorem between pure point dynamical and diffraction spectrum [5l Thm. 7], hence employing 
the diffraction measure of V calculated in [101 El] • 

3. A;-free lattice points 

The square-free integers and the visible points of the square lattice are particular cases of 
the following natural generalisation. Let A C M” be a lattice. The k-free points V = V{A, k) 
of A are then defined by 

H = A \ IJ p^A. 

p prime 

They are the points with the property that the greatest common divisor of their coordinates 
(in an arbitrary lattice basis) is not divisible by any non-trivial fcth power of an integer. 
Without restriction, we shall assume that A is unimodular, i.e. |det(T)| = 1. Moreover, we 
exclude the trivial case n = k = 1, where V consists of just the two points of A closest to 0 
on either side. On the basis of the results in [iniEi], one can then show analogous versions 
of any of the above findings. In particular, one has El Cor. 1] 

= cR) 

and the result for the diffraction measure 7 of H looks as follows. Recall that the dual or 
reciprocal lattice A* of A is 

A* := {y € M” | y • x G Z for all x G A}. 

Further, the denominator of a point i in the Q-span QA* of A* is defined as 

den(f) := gcd{m G N | mi G A*}. 

Theorem 7. |10l Thms. 3 and 5] |31l Thm. 8 ] The natural diffraction measure 7 of the 
autocorrelation of V exists. It is a positive pure point measure which is translation bounded 
and supported on the set of points in QA* with {k + l)-free denominator, so 

7 = mdi. 

ieQA* 

den(.^) (fc + l)-free 
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In particular, /(O) = (l/C(n/c))^ and when 0 7 ^ £ G Qyl* has {k + l)-free denominator den(^), 
the corresponding intensity is given by 



Again, the hull 


Xv = {t + V\teA} 

of V turns out to be just the set of admissible subsets of A, i.e. subsets A oi A with 

|Ap| <p^’^ 

for any prime p, where Ap denotes the reduction of A modulo p^A; see [311 Thm. 6 ]. The 
natural topological dynamical system (Xy,yl) has the analogous properties as in the special 
case discussed above in Theorem 01 In particular, it has positive topological entropy equal to 
the patch counting entropy of V, i.e. 

VO') = H 

by [Ml Thm. 3]. For the patch frequencies, one has the following result. 

Theorem 8. |3Tl Thms. 1 and 2] Any p-patch V of V occurs with positive frequency, which 
is given by 

.(v)= x 


J-C(AnBp(0))\7’ 


p. 


{nk 


□ 


This gives rise to a measure-theoretic dynamical system (Xy,T, p) which can be seen, as 
above, to be ergodic and metrically isomorphic to {G,A,p), where G is the compact Abelian 
group 

G = X\Ap = WA/p^A 

p p 

on which the lattice A acts via addition of l{x) = ([a)]p). As before, p is the normalised Haar 
measure on G. It follows that (Xy,A, p) has zero measure entropy; see also [311 Thm. 4]. 
Again, V turns out to be p-generic. We thus get the analogous result to Theorem [ 6 ] also in 
this more general setting. 


4. ^-FREE LATTICE POINTS 


One further generalisation step seems possible as follows. In [18] . Lemahczyk et al. studied 
the dynamical properties of .^-free systems, i.e. C {2, 3,... } consists of pairwise coprime 
integers satisfying 

1 

- < 00 
b 

b&SS 

and the hull X^ = {t + 1/ 11 G Z} is the orbit closure of the set 

y = Z \ IJ 6 Z 

b^SS 
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of SS-free numbers (integers with no factor from ^). Replacing the one-dimensional lattice 
Z C K by other unimodular lattices A C M"' in the above definitions and requiring that 
Ylb&gs W ^ arrives at SS-jree lattice points V and the associated topological dynamical 

systems (Xy, A). The A:-free lattice points from the previous section then arise from the 
particular choice ^ = {p^ \ p prime}. Since the proofs in [101 [31] do not use special properties 
of A:th powers of prime numbers except their pairwise coprimality, the above results carry over 
to the case of ,^-free lattice points with almost identical proofs. In particular, the density 
and topological (patch counting) entropy of V are given by 

dens(R) = n 

and hpc{y) = log(2) dens(R), respectively. Again, Xy contains the admissible subsets A of 
A, i.e. 

lAI < b^ 

for any b E where Aj, denotes the reduction of A modulo bA. Moreover, the diffraction 
measure 7 of R exists. It is a pure point measure that is supported on the set of points 
i E QT* with the property that the denominator den(I’) divides a finite product of distinct 
6 ’s, the intensity at such a point being given by 

1(0 = (de.,s(y) n 

\ bem 

gcd(den(£) ,b)^l 

Both the pure point nature and the intensity formula can be shown by writing 7 as a vague 
limit of diffraction measures of approximating crystallographic systems. A Weierstrafi M-test, 
as in [T 0 |, implies that one also has the stronger norm convergence, which preserves the pure 
point nature of approximating measures in the limit via |3l Thm. 8.4]. 

Further, the frequency of a p-patch P of R is positive and given by the expression 

"(p) = E (-1)'^' n 

J'c(AnBp(o))\'P 

The associated measure-theoretic dynamical system (Xy,yl, y) can be seen, as above, to 
be ergodic and metrically isomorphic to {G,A,p), where G is the compact Abelian group 

G=n =n 

be^ b£^ 

on which the lattice A acts via addition of t(x) = ([x];,) and p is the normalised Haar measure 
on G. Again, V turns out to be y-generic. Due to the similarity of the structures, we leave 
all details to the reader. 


1 - 


liV A J^)b\ 



5. fe-FREE INTEGERS IN NUMBER FIELDS 

Another possible extension is given by the number-theoretic setting that was studied by 
Cellarosi and Vinogradov in M- Let K be an algebraic number field, i.e. a finite extension 
of Q, say of degree [iL : Q] = d; see [lU |30| for background material on algebraic number 
theory. Further, consider the Dedekind domain Ok of integers in K (in particular, any ideal 
0 / a / Ok of Ok factors uniquely as a product of prime ideals) and define, for k > 2, the set 
V of k-free integers of K as the nonzero elements a E Ok with the property that, if a is not 
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Figure 3. The Minkowski embedding of Z[\/2] in 


a unit, then the prime ideal factorisation of the corresponding principal ideal 0 / (a) / Ok 
contains no feth powers of prime ideals, i.e. (a) ^ for any prime ideal p C Ok- In other 
words, one has 

V = V{OK,k) = Ok\ U P^ 

p prime ideal 

It is well known that Ok is a free Z-module of rank d and is thus isomorphic to the lattice 
as a group. In particular, there is a natural isomorphism from Ok to a lattice in namely 
the Minkowski embedding] see [I21I30] and [3l Ch. 3.4], In order to illustrate this, we prefer 
to discuss the specific real quadratic number field K = Q(\/2) with Ok = Z[\/2]. It is well 
known that K is a Galois extension and thus has precisely two field automorphisms, namely 
the identity and the non-trivial automorphism determined by \/2 i—)• — \/2- We denote the 
latter automorphism by x i-A xh One can then easily check that Ok corresponds under the 
Minkowski embedding j: iii —)• R^, given by 

X i-A (x, x'), 

to a non-unimodular lattice L in R^ with area | det£| = 2\/2 = ■\/\d^, where dK = 8 is the 
discriminant of K. In fact, since Ok = Z © Z-\/2, one has £ = Z(l, 1) © Z(-\/2, —v^) ; see 
Figure El 


Moreover, the image j(a) of any ideal 0 7 ^ o C Ok is a lattice in R^ with area 
(2) \detj{a)\=2V2iOK--a), 

where {Ok '■ a) denotes the (finite) subgroup index of a in Ok, i-e. the absolute norm N{a) 
of a. Note that the absolute norm is a totally multiplicative function on the set of non-zero 
ideals of Ok- We are thus led to consider the familiar object 

j{V) = £\ljj(p"), 

p 
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Figure 4. A central patch of the Minkowski embedding of the square-free 
integers V of Q(-\/2). 


where here and below p runs through the prime ideals of Ok] see Fig. [4] for an illustration. 
Again, the proofs in [iniEi] can be adjusted to obtain similar results as above. 


Denote by 


Ck{s) - Cq(V2)('5) - X] N(aY 

aCOK ^ ^ 
a 7 ^ 0 ideal 



—V‘ 

iV(p)V 
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the Dedekind (^-function of K, which converges for all s with Re(s) > 1. Employing Eq. ([2]), 
a similar reasoning as in the previous sections now shows that the density and topological 
(patch counting) entropy of j(V) are given by 

0(F)) = ^ ^ ^ n (l - 

and log(2) dens(j(E)), respectively. Note that the Chinese Remainder Theorem in its general 
form says that, given pairwise coprime ideals ai,..., in a ring O (ag + at = O for s / t), 
one has 01=1 a* = ai n • • • n a^. and 

r r 

O/Wat^ Wo/at. 

i=l i=\ 

Recall that the dual or reciprocal module of Ok is the fractional ideal 

^*K ■= {y ^ K I Tr^/Q(yx) G Z for all x G Ok} 

of K containing Ok, where Tr^/Q(yx) = yx + y'x' is the trace of yx. Then j{0*^) = C* and 

hence j{QO*j^) = QC*. Here, one calculates that O}^ = Zi0Z^ and thus QOk = QO}^ = K 
as well as QC* = QC. Eurther, the denominator of a point i in QT* is defined as the non-zero 
ideal 

den(£) := {x G Oio | xj-\£) G 0*k} C Ok- 

Then, the diffraction measure 7 of j{V) is pure point (for the same reason as above) and is 
supported on the set of points £ G j{QO}^) = Q£* with {k + l)-free denominator den(^) (i.e., 
either den(.£) = Ok or the unique prime ideal factorization of den(£) contains no {k + l)th 
powers), the intensity at such a point being given by 

^ n iv(p)fc -1 

den(£)Cp 

See Fig. [5] for an illustration, where the restriction of 7 to a compact region is shown. 



With a view to the general case of an algebraic number field K, the above did not make 
use of the additional fact that Z[\/2] is in fact a Euclidean domain and thus a principal ideal 
domain (in particular, a unique factorisation domain). Here, one has N{a) = |N^/Q(a)| for 
a = (a), where NK/Q{a) = aa' is the norm of a. Note that Z[\/2] is a Euclidean domain 
with respect to the norm function x i-)- |N'r/q(x)|, i.e. a + b^/2 i-A |a^ — 26^ |. Furthermore, 
the Dedekind C-fuRction of K can be written more explicitly in terms of the usual rational 
primes, i.e. 


Ck{s) 


1 


1 - 2 -^ 


n 

P=±l(8) 


1 

(1 — p “®)2 


n 

P=±3(8) 


1 

1 — p~‘^^ ’ 


see P Eq. (7)]. For instance, one has Cat(2) = ^^ 75 i see P Eq. (58)]. Hence, the intensity 
!(£} can be computed explicitly in terms of the prime elements of Ok dividing any generator 
of the (principal) ideal den(^). 

Let us finally turn to the associated topological dynamical system (Xy,Ox) — (^j(v),^), 
where Xy resp. Xj^yj are given as the translation orbit closure of V resp. j(V) with respect to 
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Figure 5. Diffraction 7 of the Minkowski embedding of F = F(Z[-v/2],2), 
with the intensities rescaled by the function x 1 —)• ^^j20 for better visibility 
of small intensities. Its lattice of periods is C* as described in the text and the 
picture shows the intensities inside the closure of one fundamental domain. 


the product topology on {0, 1}^^ resp. {0, l}"^. Again, Xy resp. '^j{v) can be characterised 
as the admissible subsets A of Ok resp. C, i.e. 

141 < N{p f 

for any prime ideal p of Ok, where Ap denotes the reduction of A modulo p^ resp. j(p^). 
Further, the frequency of a p-patch V of j{V) is positive and given by the expression 


Ffp) 


TcicnBp{o))\p p ^ 


K^u^pl A 
N{pr )' 


The associated measure-theoretic dynamical system C,u) ~ (Ky ,Ok,v) can be 

seen, as above, to be ergodic and metrically isomorphic to {G,C,fj,), where G is the compact 
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Abelian group 


0 = 114 = 11 C/jip^)c:^llOK/P^=ll{OK)p 

P P P P 

on which the lattice L (resp. the group Ok) acts via addition of i{j{x)) = ([j(x)]p) (resp. 
i{x) = ([x]p)), where x G Ok, and /i is the normalised Haar measure on G. As above, V turns 
out to be p-generic. 

Since none of the above uses special properties of the quadratic field Q(\/2), similar results 
hold for the general case of an arbitrary algebraic number field K. Moreover, even the 
extension to ,^^-free integers in K, i.e. 


V = Vo^ = OK\^b, 

b&SS 


where is a set of pairwise coprime ideals b C Ok satisfying 


E 


1 

W) 


< oo. 


should be possible. 


Appendix A. Ergodicity of the patch frequency measure 

Below, we shall only treat the paradigmatic case V = V ^2 from Section [2l For a p-patch 
P G A(p) of V (i.e. P C Bp(0) H with P G A), denote by Bp the set of elements of 
Xv' = A whose /9-patch at 0 contains P. One readily checks that the sets of type Bp form a 
semi-algebra that also generates the Borel c-algebra on Xy- In fact, one has 

Cp = Bp\ Bq, 

SeA(p) 

V<ZQ 


and 


(Al) Bp= [j Cq. 

QeA(p) 

VcQ 


Corollary 1. For any p-patch P ofV, one has 


V 


{Bp) = 11(1 
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Proof. Let V he a p-patch of V. In this proof, we indicate summation variables by a dot 
under the symbol. By Theorem [5] and the definition of i', (HH) implies that 




E ^(^2) 


Q&A{p) 

VCQ 


E E 1 


Q&A{p) Q(iTcJ?r\Bp{o) 
VCQ 


\^p 

p2 


E 


VcQcTcZ^r\Bp{0) 


\^p 

p2 


VcQcVc'z'^nBpio) p 

T+V 




n 1 


Vv 


p 


,2 / ’ 


since, for fixed C n Bp{0) with V C B, one indeed has 


E 1 - 

VCQCT p ^ 





VCQCT 



where the last equality follows from the binomial theorem since J-\V 7 ^ 0 . 


□ 


For a natural number m, finite subsets V and Q of 1? and S C Vm-, we set 

~ ( ri ~ '^) ^ ( U ~ '^) ’ 

SCS SCVm\S 

i.e. the set of elements of (Z^)m that lie in Qm — s precisely for those s G Vm with s G 

S C Vm, in particular Q^-p = (Z^)m \ (Qm — Vm)- With q'^p = |Q^-p|, one then has 

q^P = m? — \ Qm — Vm\ and, since the difference set Qm — Vm is the disjoint union of the 
various Q™-p, where one has 0 / S' C Vm, 

(A2) E 

SCPm 

Note that the following two lemmas also hold for any finite subsets V and Q of an arbitrary 
finite group G instead of G = {7?)m- 
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Lemma 2. For any finite subsets V and Q of 7? and any natural number m, one has 

Y. \S\(isF = \'^rn\\Qm\. 

SdVm 

Proof. We use induction on \Vm\- For the induction basis \Vm\ = 0, the assertion is trivially 
true. For the induction step, consider V with \Vm\ > 0 and fix an element, say *, of Vm- It 
follows that 


E = E E (i^i + i)95uw,p 

SCVrr, SCPmM*} SCVm\{*} 

= E \smv+<isu{.},v) + E 

SCP„r\{*} SCPm.\{*} 

= E + |Qm|, 


Since 


E «?□(.),P = E 

SCPm\{*} SCPm\{*} 


\ (^ y Qm- S 
seP„,\{Su{*}) 


( {Qm {*}) Fl ^ y Qm 

V se 5 

{Qm - {*}) n f Pi Qm - s) \ ( y Qm- sj 


= E 

5 cP™\{*} I \ ses 

— I Qm ~ {*} I 
= |Qm|- 

Furthermore, for S C Vm \ {*}> one has 

Qsum,v u Qtv 

{Qm - {*}) n f P Qm - s) \ ( y Qm- sj 


se{Vm\{*})\S 


se5 




U 


((n2™-«)\(fi™-w))n((n©»>-»)\( u fl”-* 


(nc- 


— s 


ses 


\y y Qm- S 
S&iPm\{*})\S 


— (O™ 

- ^S,P\{*} 

wherefore, by the induction hypothesis, 

E IS|(<iSp + ' 7 ?u{.>,p)= E 

SCPruM*} SCPm\{*} 

This completes the proof. 


□ 
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Lemma 3. For any finite subsets V and Q ofl? and square-free d, one has 

n ^ ^ 5 ,^) = W_{P^\"Pp\F P^\Qp\ ~ \"Pp\\Qv'\) ■ 

{vp)p\d p\d ScVp p\d 

0<Up<\Vp\ \S\ = \Vp\-Up 

Proof. This is proved by induction on the number u}{d) of prime factors of d. For the induction 
basis u}{d) = 1, say d = p, note that by (IA2I1 and Lemma [2] one indeed has 

((p'p + |2pl) ^ ^S,p) ~ (|2pl ^ 9s,-p) + (I'^pl ~ I^^D^g.-p 

0<i'p<\'Pp\ SCPp ScVp ScVp 

\S\ = \Pp\-Up 

— |Qp|p P I'Pplp “I'^pllQpl- 


For the induction step, consider d with uj{d) = r +1, where r > 1, say d = pi.. .prPr+i. Then 


n((^p+i 2 pi) X/ ^s,p 

Mp\d p\d ScVp 

0<Up<\Vp\ ISHVpl-i'p 

Ppr+l I 


(* + lQpr+ii) ^ <?5y-p^ ((p'p + i 2 pi) ^ Qs, 


*—0 {p'p)p\d 

0<t'p<|'Pp| 

IP’pr+l I 

^ (* + I Qpr+1 I 

i = 0 


\S\ = \Vp^+f-i 




ScVp 

\S\ = \Vp\-Up 


S^T^Pr+l 

\sHVp^+f-i 


^ ((p'p+iQpi) X] 

Mpl^pl^ 


'Pr +1 Pr+l 

0 <^'p<i'Pp| 


ScVp 

\S\ = \Vp\-yp 


which is, by the induction hypothesis. 


Ppr+l I 

y~! (^ + 12; 

j=0 


■Pr+1 I 


E <i%v n {p>p\+p^\Qp\-\'Pp\\Qp\) 


Scv, 


Pr+1 


\S\ = \Vp^+d-i 


P\ 


Pr+l 


I’P’pr+ll 


n (p^l^pl+p^l2pl - l^pllQpl) E E (^ +I2p.+il)9: 


,Pr+l 

s.-p ■ 


Pl 


j =0 Scv, 


Pr + 1 


Pr+1 


|S| = |Pp,+i|-i 


It thus remains to show that 

\Ppr+l I 

El El (^ + |2pr+ll)95,p = Pr+ll’Ppr+ll +Pr+l|2pT-+ll “ l"^Pr+ll|2pr-+ll) 

i =0 SCPp ^+1 

|5| = |Pp,+i|-* 


ERGODIC PROPERTIES OF VISIBLE LATTICE POINTS 


21 


which is clear since, by (|A2p and Lemma [2] again, one has 


Tpr + l 

1 





E 



\Q-Pr+l\)^sfp 



i=0 






|5| = 

[Ppr+l\ ® 




\Ppr+l 

1 

\Ppr+l 1 



= 

E 

E 


E 

2pr+l 


j=0 

SC.Ppr+1 

i=0 





|S| = |Pp,+ i|-l 


|S|=IW.+il- 

—i 

= 

E 

(l"^Pr+l 1 ~ 

\^\)Qr,P^ + Qpr+l ^ 

Pr+l 

^S,P 




~ \'^Pr+l\Pr+l I I 2pr+l I + I 2pr+l |Pr+l- 


This completes the proof. 


□ 


The proof of the main result below needs the well-known estimate 
(A3) \Bp{x)f\I?\=T:p^+ 0{p) + 0{l), 

approximating the number of points of in a large ball Bp{x) (the last error term being 
required only when < 1). This is obtained by dividing Bp{x) into fundamental regions 
for 1?, each of volume 1 and containing one point of 1?, with the error terms arising from 
fundamental regions that overlap the boundary of Bp(x). A more precise version is given 
as [ini Prop. 1]. 

Theorem 9. The measure v is ergodic. 


Proof. We have to show that 


v{B-p)v{Bq) 


1 


lim 
R-^oo vri?^ 


^ v{x + Bpr\ Bq) 


for any finite admissible subsets V and Q of Z^; cf. [411 Thm. 1.17(i)]. By Corollary [H the 
latter is equivalent to 


n 




1 


lim 
R-^oo vr/t^ 


E 

3 ;SS_R( 0 )nZ 2 


p(x -I- Bp n Bq). 


Let p denote the Mobius function for the rest of this proof. The left-hand side is equal to 


V V p p p 


p V \p p p 

d square-free p\d ^ 

E ^n(iAi+iQ,i- 

d square-free p\d ^ 


\'Pp\\Qp\ \ 

J 

p 2 j ■ 
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For a fixed R > 0, the right-hand side is 

1 ^ VI + 

ttR^ lir p2 

— y 

-d2 


ttR? 

1 

-kR? 

1 




\{QVJx + V)^ 

p2 


xeBji(0)nZ^ d square-free p\d 

d square-free x&Bii(0)OJ? p\d 

L E E 




E n 

d square-free Wp)p\d x&Bii{0)fM? p\d 

\Qp\<Vp<\Qp\+\Vj,\ 'ip\d-.\{QVJx+V)p\=Vp 

w E M(^) E n5 E 




ttR? 


1 


pz 

d square-free (i^pipld pM xe-B_R(0)nZ2 

|Spl<i^P<|Sp|+|Pp| Vp|d:|(SUx+7’)p|=P'p 


— E 

-D2 


^(d) 


E E 


7ri?2 d? 

d square-free {i^p)p\d P\d x&Bii{0)r\I? 

\Qp\<Vp<\Qp\+\Pp\ 'ip\d-.\{QVJx+V)p\=Vp 


1 . 


By |10l Prop. 1] and Proposition [2] with 1? replaced by the lattices pi?, where p \ d, one 
obtains the estimate 

,2x1-1/2\ \ 

I + 0(1) I \{xp € {l‘^)p : |(Q U X + V)p\ = i'p}\ 

/ ) p\d 


'kBt 

~?P 


+ o 


y 


for the inner sum. Substituting this in the above expression and letting R tend to inhnity, 
one obtains 


E ^ E Yiyiixp^i^p ■ \{Q^x+v)p\ = up}\) 

d square-free (*^p)p|d p\d 

\Qp\<iyp<\Qp\ + \Pp\ 

^ E ^^4^ E n ((^P ■*“ l^pl^ K^P ^ O a:-I-P)p| = t'p-I-|Qp|}|) 

d square-free {'^p)p\d P\d 

0<i/p<|'Pp| 

and the inner product can be rewritten as 

{{vp + IQpl) |{xp E {1^)p : Xp E Q^g-p for S C Vp with l^l = |Pp| - z/p}|^ 
p\d 

~ n((^p + |2pl) E 

p\d ScPp 

\s\=\rp\-up 

Thus, in order to prove the claim, it suffices to show, for square-free d, the identity 

E n(K + l2pl) E = IIEl^pl+p'IQpl-l^pllQpl), 

(i^plpld p\d ScVp p\d 

0<Vp<\Vp\ \S\ = \Pp\-Vp 
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which is just the content of Lemma [3l □ 

Appendix B. Isomorphism between (Xy,Z^,p) and {Y\p{Z‘^)p,Z‘^, ^) 

Let Ai be the Borel subset of Xy = A consisting of the elements X G Xy that satisfy 

\Xp\=p^-l 

for any prime p, i.e. X misses exactly one coset of in Z^. There is a natural Borel map 

9: Ai ^T, 

given by A I—>• {[yp]p), where [yp\p is nniquely determined by [yp\p 0 Xp. Note that 9 fails to 
be continnons or injective. Clearly, Ai is Z^-invariant and 9 intertwines the Z^-actions. Note 
also that, for A G Ai, one clearly has 

(Bl) Xcpi9iX)). 

Lemma 4. One has IL g Ai. 

Proof. Fix a prime nnmber p and choose a set Aoip^ — 1 elements of Z^ snch that \Ap\ = —1. 

By the Chinese Remainder Theorem (Prop. [2]), we may assnme that Apt = {(0,0)p'} for the 
finitely many primes p' < p. Then, A G A and, by Theorem [SJ there is a translation f G Z^ 
such that f + A C V. Since 

- 1 = |(t + A)p| < \Vp\ < - 1, 

the assertion follows. □ 

Lemma 5. One has p(Ai) = 1. 

Proof. The ergodicity of the fnll measnre v (Thms. [5] and [ 6 Kb)) implies that y-almost every 
A G Xy has a dense Z^-orbit; compare im Thm. 1.7]. It follows from Lemma 0] that 

{A G Xy I A has a dense Z^-orbit} C Ai. 

This inclusion is dne to the fact that, if A has a dense orbit, then in particnlar P G Ai is 
an element of the orbit closure of A. Since, for any prime p, representatives of the p^ — 1 
different elements of Vp in V can be chosen within a finite distance from the origin, it is clear 
from the definition of the topology on Xy that there is a translation t G Z^ snch that t + A 
contains these representatives. Thns 

p^ - 1 > |Ap| = |(t + A)p| > |Pp| =p^ - 1 

and the assertion follows. □ 

Hence the push-forward measure of to T by the map 0 is a Z^-invariant probability 
measure and thus is the normalised Haar measnre p. Set 

Ti := (/?"^(Ai), 

which is a Z^-invariant Borel set with <^(Ti) C Ai. In particular, this shows that 9o(p\f.^ = idxi 
and thus the restriction of p to Ti and the restriction of 9 to (/^(Tri) are injective. 

Lemma 6. One has 0(Ai) = Ti. 
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Proof. It suffices to prove the inclusion 9{Ai) C Ti, since this immediately yields the assertion 
due to Ti = 0((/?(Ti)) C 9{Ai). So let us assume the existence of an X E Ai with 9{X) 0 Ti, 
i.e. ip{9{X)) 0 Ai. Then there is a prime number p such that 

IM0(x)))p| <p2_i. 

Using (IBip . this implies that also \Xp\ < — 1, a contradiction. □ 

In particular, this shows that 0“^(Ti) = Ai and thus 

/r(Ti) = v{Ai) = 1. 

It follows that 0 is a factor map from to , p). 

In order to see that 9 is in fact an isomorphism, let us first note that the Borel map 
(/j: T —)• Xy is measure-preserving since, for any p-patch V, one has 

= p({([yp]p) G TI [yp]p 0 Vp for all p}) 



= v{Bp) 


by Corollary [TJ Next, consider the subset A^ of elements X E Ai that are maximal elements 
of A with respect to inclusion, i.e. 

VU E A : (X C y ^ X = y). 

Clearly, V and every translate of V are elements of A^; see Lemma [H Using Lemma El one 
can verify that A^ contains precisely the elements X E Ai with 

(B2) X = ^(0(X)). 

Employing (|B2p . one further verifies that 

AI = Ain^(T). 

Since <p(T) can be seen to be a Borel set, Al| is thus a Borel set with measure 
u{Al) = i/(Ai n (p(T)) = i/(Ai) + u{p{T)) - y(Ai U <p(T)) = I. 

Setting 

T^ ■.= p-\Al), 

the restrictions p: T^l —)• Al| and 9: A^ —T^ are well-defined (X E Al| => (p(0(X)) = X E 
All) and can now be shown to be bijective and inverses to each other. Hence 9 and p are 
isomorphisms. 
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